THE (UNSTABLE) THRESHOLD OF BLACK HOLE FORMATION 



M.W. CHOPTUIK 

Center for Relativity 
Dept. of Physics 

The University of Texas at Austin 
Austin TX, 78712-1081 



1. Introduction 

In recent years it has become apparent that intriguing phenomenology ex- 
ists at the threshold of black hole formation in a large class of general 
relativistic collapse models. This phenomenology, which includes scaling, 
self-similarity and universality, is largely analogous to statistical mechan- 
ical critical behaviour, a fact which was first noted empirically, and sub- 
sequently clarified by perturbative calculations which borrow on ideas and 
techniques from dynamical systems theory and renormalization group the- 
ory. This contribution, which closely parallels my talk at the conference, 
consists of an overview of the considerable "zoo" ' of critical solutions which 
have been discovered thus far, along with a brief discussion of how we cur- 
rently understand the nature of these solutions from the point of view of 
perturbation theory. The reader who wishes additional details concerning 
the subject is referred to Gundlach's excellent recent review [1]. Earlier syn- 
opses by Evans [2] and Eardley [3] may also be of interest. Finally, those 
readers interested in the relationship of black-hole critical phenomena to 
questions of cosmic censorship — a topic which is not discussed below — can 
consult Wald's recent discussion of the status of cosmic censorship [4]. 



2. Critical Behaviour in Massless Scalar Collapse 

I begin with a fairly detailed description of critical behaviour in the self- 
gravitating dynamics of a spherically-symmetric, massless scalar field. Apart 
from historical reasons, I do so, not only because it is arguably the simplest 
"realistic" model of gravitational collapse, but, more importantly, because 
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most of the features we currently associate with critical collapse can be 
seen in the model. 



2.1. EQUATIONS OF MOTION 

Adopting the usual spherical coordinates (t,r,6,(p), and geometric units, 
G = c = 1, the spherically-symmetric spacetime metric can be written 

ds 1 = -a 2 (r, t) dt 2 + a 2 (r, t) dr 2 + r 2 (d9 2 + sin 2 9 dip 2 ) • (1) 

The coordinate system I have thus adopted is a natural generalization of 
Schwarzschild coordinates — in numerical relativity parlance it is the polar- 
areal (or polar-radial) system. Note that the radial coordinate, r, measures 
proper surface area and thus has an immediate geometric interpretation. 
The time coordinate, t, on the other hand, has no particular physical sig- 
nificance. However, from the point of view of critical phenomena, there is a 
preferred labeling (reparametrization) of the t = constant surfaces, namely 
the one given by the proper time, To(t), of an observer at rest at r = 0: 

T (t) = f t a(0,t)dt. (2) 

Jo 

Defining auxiliary scalar field variables, $ and II: 

*(r,f) s ^(r,t), (3) 

n(M) , ^fM), (4) 

a sufficient set of equations for the EMKG (Einstein-massless-Klein-Gordon) 
model is 
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Here (7) is the slicing condition which constrains the lapse function, a, at 
all instants of time, and (8) is the Hamiltonian constraint which similarly 
constrains the radial metric function a. 
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It is useful to discuss the dynamics of the scalar field in terms of further 
auxiliary variables, X and Y, defined by 

X(r,t) = s/2^-$ = V2tF-^, (9) 
a a or 

Y(r,t) = y/2^-U = y/^-^. (10) 
a a ot 

Note that the equations of motion (5)-(8) are invariant under the trivial 
rescalings r — ► kr, t — > fci, for arbitrary k > 0, and that X and Y are 
form- invariant under such transformations. Is is also convenient to intro- 
duce the mass aspect function, m(r,t), defined in analogy with the usual 
Schwarzschild form of the static spherically-symmetric metric: 



In terms of the X and Y variables, dm/dr = X 2 + Y 2 , and the total mass 
(ADM mass), Madm, of the spacetime is given by 

roo A™ roc 

M ADM = / ^dr= / X 2 + Y 2 dr. (12) 
Jo dr Jo 

Another useful relationship expresses the spacetime curvature scalar, R, as 
a function of X and Y: 

R = -8itT = 8^V>V^ = ^ (X 2 - y 2 ) . (13) 

where T is the trace of the scalar field stress-energy tensor. 

As is well known, polar-areal coordinates cannot cross apparent hori- 
zons, and thus, for the most part, cannot penetrate event horizons. However, 
black hole formation is clearly signaled in a calculation (see Figure 1) by 
(among other things) 

- -1, (14) 

r Rbh 

at some radius r = i?BH from which the mass, Mbh = -Rbh/2, of the final 
black hole can be quite accurately estimated. 



2.2. COMPETITION AND THE THRESHOLD OF BH FORMATION 

At a heuristic level, the existence of critical behaviour in the EMKG sys- 
tem and other models is a direct result of competition in the dynamics. The 
nature of this competition can be seen by addressing the question: "given 
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Figure 1. Behaviour of the mass aspect function, m(r, i), for cases of complete dispersal 
(top), and black hole formation (bottom), in massless scalar collapse. The origin, (0,0), 
in both plots is at the extreme right. Initial data in both cases is a single ingoing Gaussian 
profile — 4>{r, 0) = p g(r) — the calculations differ only in the choice of the overall amplitude 
factor, p, which is sub-critical (p < p*) in the first instance, and super-critical (p > p*) 
in the second. Note that flat regions in the plots are vacuum. In the sub-critical case, the 
scalar field implodes through the origin, then completely disperses leaving (essentially) 
flat spacetime in its wake. In the super-critical case, the scalar field again implodes 
through the origin, but now forms a black hole containing roughly half of the total mass 
of the spacetime. Dynamically, the geometry within the dispersing (outgoing) scalar field 
settles down to an exterior-Schwarzschild solution on a time scale set by the size of the 
hole, T ~ .Rbh = 2M BH . 
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Figure 2. Schematic illustration of black-hole-threshold possibilities. The top panel 
represents a first-order (Type I) transition, where the "order parameter" (the black hole 
mass) exhibits a gap at threshold, while the bottom panel shows a second-order (Type 
II) transition, where the black hole mass is infinitesimal at the critical point. 

generic initial data representing an imploding pulse (shell) of scalar radia- 
tion, where can the energy in the system end up at late times?" Roughly 
speaking, the kinetic energy of the massless field wants to disperse the 
field to infinity, whereas the gravitational potential energy (entirely self- 
induced), if sufficiently dominant during the collapse, will result in the 
trapping of some amount of the mass-energy of the system in a black hole. 
Indeed, the fact that, for generic initial data, there are only these two qual- 
itatively distinct end-states in the model has been rigorously established 
by Christodoulou [5, 4]. However, empirical evidence (i.e. direct solution 
of the equations of motion) rapidly leads one to the same conclusion. The 
key point is that the dynamical competition can be controlled by tuning a 
parameter in the initial conditions: it is an easy matter to set up families of 
initial data, {<J?(r, 0;p), LT(r, 0;p)} such that if the parameter, p, is less than 
some critical (threshold) value, p*, the scalar field completely disperses, 
while if p > p* a black hole forms (Figure 1). 
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Figure 3. Near-critical evolution of the auxiliary scalar field variable, X (see equa- 
tion (9)), in massless scalar collapse. Evolution proceeds left-to-right, top-to-bottom. In 
each frame, X is plotted versus the logarithmic radial coordinate, p = lnr, with dis- 
played ordinate values ranging from -0.37 to 0.33. Constant increments of logarithmic 
time, r = ln(r ( J — To), separate each successive snapshot. The discrete self-similarity of 
the critical solution can be clearly seen by comparing the trailing (left-most) edges of 
pairs of waveforms separated by two rows in the plot. Each such pair is very nearly in the 
same "phase" of critical evolution, but the dynamics of the later frame (more negative r) 
is occurring on a scale some exp(A) ~ 30 times smaller than that of the earlier snapshot. 
Note that the leading parts of the waveforms (the roughly sinusoidal oscillations) are 
almost purely-outgoing, and do not, per se, constitute part of the critical dynamics. The 
fact that these oscillations appear to be "frozen" , rather than propagating, is a result of 
the fact that the time between successive frames is exponentially decreasing. At any scale 
(position along the horizontal axis) of the critical evolution, there are three basic possi- 
bilities: (1) if the evolution is sub-critical (at that scale), the scalar field will completely 
disperse, (2) if the evolution is super-critical, a black hole with a size set by the current 
scale of the critical dynamics will form, and (3) if the evolution is precisely-critical, the 
strong field evolution will continue to a smaller scale. This behaviour provides convincing 
evidence for a Type II transition in the model. Plots of the other scalar field variable, 
Y, as well as X 2 + Y 2 — dm/dr and X 2 — Y 2 — 4R/r 2 show analogous oscillatory be- 
haviour between fixed limits. Since the strong-field regime is characterized by r — > for 
a precisely critical solution, this last fact shows immediately that the spacetime scalar 
curvature, R, diverges in the critical limit. 
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The existence of interpolating families and black hole thresholds in the 
EMKG model was well-established in early studies of the system [6, 7, 8]. 
The discovery of critical behaviour in the model, however, came later, and 
was a direct result of a question Christodoulou posed in 1987 [9]: "will 
black hole formation turn on at finite or infinitesimal mass for a generic 
interpolating family at threshold?" (see Figure 2). Viewing the black-hole 
mass as an order parameter, I refer to these two possibilities as Type I and 
Type II transitions, respectively, in analogy with first and second order 
phase transitions in statistical mechanical systems. 

2.3. TYPE II BEHAVIOUR IN THE EMKG MODEL 

Detailed phenomenological studies [10, 11, 12] using many initial data fami- 
lies clearly demonstrated that the strong-field (i.e. critical) dynamics in the 
EMKG model is characterized by an essentially unique solution of the equa- 
tions of motion, which I denote schematically as Z*. Like the Schwarzschild 
solution, Z*, is only determined upto an overall (length/time) scale, but 
unlike Schwarzschild, the critical solution is obviously (by construction!) 
unstable — the slightest perturbation will result in either complete dispersal 
or black hole formation. 

A key feature of the critical solution (and of all currently known Type II 
solutions) is self- similarity or scale invariance — loosely speaking, the criti- 
cal solution has a scale (homothetic) Killing vector. (Correspondingly, the 
known Type I critical solutions discussed in section 6 are characterized by 
timelike Killing vectors). Empirically, this means as one tunes closer and 
closer to a critical point one sees, in a single evolution, the same strong-field 
dynamics playing out on a wider and wider range of spatial scales, always in 
a (shrinking) neighborhood of r = 0. In the scalar field case (see Figure 3), 
the dynamics at any particular "scale epoch" is decidely non-trivial — i.e. 
the self-similarity is discrete, rather than continuous. This novel "echoing" 
feature of the solution, whereby the dynamics repeats on scales related by 
a factor of e A , A = 3.44..., is one of the most intriguing aspects of the 
solution, and its origin remains somewhat of a mystery. We can express the 
discrete self-similarity is slightly more mathematical terms by noting that 
in a precisely critical evolution, the strong-field evolution "accumulates" at 
a singular event with (r, To) coordinates (0, Tg). Then defining logarithmic 
coordinates, p = lnr, r = hi(Tg — Tg), we have 

Z* (p ± nA, r ± nA) ~ Z* (p, r) ra = 0, 1,2, ■ ■ ■ . (15) 

As explained in more detail in Figure 3, the nature of the critical dy- 
namics (which is now well-understood, albeit in a non-rigorous fashion, 
in perturbation theory — see section 5), makes it clear that tuning of ini- 
tial data provides the mechanism for making arbitrarily small black holes, 
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Figure 4- Typical evidence for mass scaling in the collapse of a massless scalar field. 
The data displayed here was generated from a one-parameter family of initially-ingoing 
Gaussian pulses of scalar field in which the overall amplitude of the pulse was varied. The 
mass-scaling exponent, 7ls ~ 0.376, was determined from a least-squares fit. The inset 
clearly shows the "second order" (Type II) nature of the transition. It is worth noting 
that, particularly for compact pulse shapes such as the Gaussians used here, the mass 
scaling persists well out of the asymptotic regime, p — > p* . 

which in itself establishes that Type II transitions occur generically in the 
EMKG model. Direct measurement of the black hole masses in the super- 
critical regime provides additional evidence. One finds that, for a generic 
interpolating family S[p], and as p — > p*, the black hole masses are well- fit 
by a scaling law (see Figure 4): 



M BH = c f \p-p* 



(16) 



THE (UNSTABLE) THRESHOLD OF BLACK HOLE FORMATION 9 



where c/ is a family-dependent constant, but 7 = 0.37... is universal (i.e. 
family- independent) . 

As also noted in Figure 3, the nature of the critical solution immedi- 
ately implies that the scalar curvature (see equation (13)) grows without 
bound near r = in a precisely critical evolution. It is also clear from the 
simulations that these regions of arbitrarily high curvature are visible by 
observers at infinity. Indeed, the trailing edge of the scalar field pulse it- 
self effectively does the job of tracking outwards propagating null geodesies 
emitted from the critical region. Again, the point is that, contrary to some 
expectations [13], without the formation of a horizon, one cannot indefi- 
nitely localize the mass-energy of the scalar field. Finally, I note that all 
of the above features of critical spherically-symmetric scalar collapse have 
been reproduced — using a variety of theoretical approaches and numerical 
techniques — by several other groups [13, 14, 15, 16]. 

3. Other Early Results in Critical Collapse 

3.1. GRAVITATIONAL WAVE COLLAPSE 

Shortly after the discovery of critical behaviour in the EMKG model, Abra- 
hams and Evans [17] presented exciting results demonstrating that very 
similar effects occurred in the purely gravitational (i.e. vacuum) collapse of 
axisymmetric gravitational waves. Again, at least in principle, it is straight- 
forward to construct interpolating families in this case, but, computation- 
ally, this problem is much more difficult to treat than the spherically sym- 
metric EMKG system described above. Nonetheless, Abrahams and Evans 
were able to produce convincing evidence for (1) a discretely self-similar 
threshold solution, this time with A ~ 0.6, and (2) mass-scaling in the 
super-critical regime, with 7 « 0.37. The work in [17] used a single inter- 
polating family; a follow-up paper [18] reported the observation of similar 
results from a second family, thus providing evidence for the universality 
of their critical solution. 

3.2. RADIATION FLUID COLLAPSE 

A third and very important example of critical behaviour came with the 
work by Evans and Coleman [19], who studied the spherically symmetric 
collapse of a perfect fluid with the simple equation of state 

P = \p, (17) 

where P and p are the fluid's pressure and density respectively. Early on, 
Evans had realized that the self-similar nature of the Type II transitions 
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which were being observed could be used to great advantage in understand- 
ing what was happening in critical collapse. Familiar with the existence of 
continuously self-similar solutions in relativistic fluid flow, he argued that 
that the phase transition in the case of radiation fluid collapse should be 
characterized by continuous self-similarity. Proceeding from the ansatz of 
self-similarity, Evans and Coleman constructed a precisely self-similar so- 
lution and verified that it was identical to the critical solution generated 
from full dynamical evolution (again using interpolating families). Once 
more, convincing evidence for mass-scaling in the super-critical regime, with 
7 0.36, was found. Importantly, Evans also suggested that an investiga- 
tion of the perturbative mode structure of the critical solution could provide 
the basis for the computation of the mass-scaling exponent. Finally, the ap- 
parent numerical equality of 7 for the EMKG, vacuum axisymmetric and 
radiation fluid cases, led to a provocative, but short-lived conjecture, that 
there might be "true" universality of 7 across various collapse models [19]. 

4. Self Similarity in Critical Collapse 

As Evans anticipated, the self-similar nature of Type II critical solutions has 
proven crucial to our current understanding of black hole critical phenom- 
ena. Here I will only briefly sketch some key ideas concerning self-similarity 
in this context — the interested reader can consult Gundlach's review [1] 
(and references contained therein) for a much more thorough presentation. 

Restricting attention to spherical symmetry, we choose "geometric" co- 
ordinates, R and T, which are adapted to any particular critical solution. 
In the context of polar-radial gauge discussed above, obvious choices are 
R = r (areal radius) and T = (Tq — To), so that the (0,0) is the central 
singularity of the precisely critical solution, and time advances to the past 
of the singularity. The natural similarity variable (coordinate) for a type II 
critical solution is then 



and one can use either T or R, or more conveniently, r = InT or p = lni? 
as a second coordinate. Evolution along lines of constant ( (see Figure 5) 
then represents evolution in scale. Continuously self-similar (CSS) solutions 
are characterized by 



where g denotes any dynamical variable which exhibits scaling. Discretely 
self-similar (DSS) solutions, on the other hand, satisfy 




(18) 



g((,T)=g{(,T f ) 



t , t' arbitrary, 



(19) 



5(C,t) = g((,T±nA) 



n = 0,1,2, ••• , 



(20) 
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p = In R 



Figure 5. Schematic illustration of continuous (CSS) and discrete (DSS) self-similarity 
in critical collapse. The similarity variable £ = R/T is constant along dashed lines. In the 
case of CSS evolution (left part of the plot), a single critical profile propagates to the left 
(i.e. to smaller spatio-temporal scales), while in the DSS case, a non-trivial, but periodic, 
sequence of profiles, is generated. The periodicity in log-space, A, is a model-dependent 
(but initial-data- independent) scaling exponent. 

where A is the model-dependent "echoing exponent" . It is important to note 
that, in both cases, the self-similar solution is not a complete spacetime, 
but (in R, T coordinates) a wedge-shaped region which, one expects, can 
be analytically continued (see [20, 21]) to produce a complete spacetime. 

5. Perturbation Theory: Mode Structure of Critical Solutions 

Again, as Evans suggested [19], the application of perturbation theory to 
black-hole threshold solutions has proven to be a powerful technique for 
studying and explaining many of the key features of critical collapse. Briefly, 
for any given collapse model (which includes specification of the matter 
content, coupling to gravity and other fields, symmetries, etc.), all of the 
current evidence suggests that one will find one or more isolated (in "initial- 



12 



M.W. CHOPTUIK 



data space") critical solutions, g*. Operationally, g* can be constructed 
either (1) indirectly, using the partial differential equations of motion and 
interpolating families of initial data, or (2) directly, starting from an ansatz 
which reflects the particular type of self-similarity exhibited by the critical 
solution. 

Schematically (see [1] and [22] for more details), one then considers per- 
turbations about the background, g*, and looks for eigenmodes, /j(C) with 
corresponding eigenvalues, Aj, with the assumption that (small) departures, 
5g(C,T), from the critical solution can be then be expressed as 

Sg(C, r) = g(C, t) - g*(C, r) = £ d exp (X t r) f t (C) , (21) 

i 

where the C, are coefficients. 

In terms of explaining the mass-scaling and universality of the Type 
II transitions discussed above, the crucial observation was made by Koike, 
Hara and Adachi [23], who noted that the "sharp" nature of the transition 
(as well as the universality) strongly suggested that of all the modes, only 
one, /*, has an associated eigenvalue, A* with 

ReA*>0. (22) 

That is, there is only one growing mode associated with the critical solution — 
all other modes generically die off as the critical solution is approached. It 
is then a straightforward matter of linearization of (21) and dimensional 
analysis to show that the black hole mass will satisfy a scaling law 

M BH «|p-p1 1/RcA * , (23) 

so that 

l = lik- < 24 > 

This argument provides an immediate and intuitive explanation for the 
observed universality (initial-data independence) of the Type II critical 
solutions — in tuning the family parameter, p, to the threshold value p*, 
one is effectively tuning out the unstable mode in the initial conditions, 
so that as p — > p*, the strong- field dynamics is well-described by the pre- 
cisely critical solution over more and more decades of scale. Moreover, the 
eventual departure from g* in any specific computation (sub-critical or 
super-critical) is well-approximated by 

fy(C,T)~exp(A*r)/*(C), (25) 



independently of any specifics of the initial data. Thus, although Type II 
critical solutions are clearly unstable, they are, in some sense, the least 
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unstable solutions possible, by virtue of the fact that they possess only a 
single growing mode. 

Koike et al tested their ideas in the context of the radiation fluid critical 
solution, found strong evidence that there was only one growing mode in 
perturbation theory, and from the eigenvalue of the mode, computed a 
scaling exponent, 7 = 0.3558 ... in excellent agreement with the Evans and 
Coleman simulations. Concurrent with [23], Maison [24] presented similar 
calculations which were, moreover, generalized to the case 

P = kp 0.01 < k < 0.888. (26) 

His specific results for k = 1/3 were in precise agreement with Koike et 
al, and, moreover, his general results clearly showed that the mass-scaling 
exponent was dependent on k — ranging from 7 = 0.114 for k = 0.01 to 
7 = 0.852 for k = 0.888. This, of course, provided strong evidence for the 
model- dependence (i.e. "non-universality") of 7. 

The picture of Type II critical solutions as "intermediate attractors" [25] 
with 1-dimensional unstable manifolds has now been validated for many 
spherically symmetric models, including axion/dilaton collapse [26, 27] (see 
also [28] for dynamical evolutions in a generalized setting), the EMKG 
model [20, 29], the SU(2) Einstein- Yang- Mills model briefly discussed be- 
low [30], and scalar electrodynamics [31]. In the last case, Gundlach and 
Martin-Garcia also predicted an additional scaling law for the (scalar) 
charge of black holes, which was very quickly observed in independent work 
by Hod and Piran [16]. 

In addition to the critical solutions with one unstable mode described 
above, other self-similar, strong field solutions have been constructed and 
investigated (see for example, the work by Hirschmann and Eardley [21, 
32] ) . It seems reasonable to assume that such solutions also sit at the black 
hole threshold, but they generically have two or more growing modes, so 
will not be intermediate attractors for generic 1-parameter interpolating 
families. However, there is evidence [33] that at least some of these solutions 
may be constructed via carefully constructed families of initial data. 

Gundlach's perturbative analysis of the EMKG system [20] also lead him 
to predict the existence of a universal "wiggle" in the mass-scaling law (16). 
Hod and Piran presented similar arguments and numerical evidence for the 
oscillation in [15]. It is probable that this wiggle is visible in Figure 4, but 
the systematic errors in the finite-difference computations which produced 
those results preclude a definitive statement. 

Finally in interesting recent work, Garfinkle [34] considered the separa- 
tion of Einstein's equations into pieces governing (1) the dynamics of the 
overall scale of a self-gravitating system, and (2) the dynamics of the "scale 
invariant" part of the metric. As part of his study he constructed a model 
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Figure 6. Evolution of the Yang- Mills potential, W(r), in a near-critical evolution at 
the Type I transition in Einstein- Yang-Mills collapse. Solid lines show the dynamical evo- 
lution of W(r), while the dashed lines superimpose the static, n = 1, Bartnik-Mckinnon 
solution. The initial data is an incoming "kink" which implodes to the center (frames 
1-3), sheds off Yang-Mills radiation and asymptotes to the static solution (frames 4-9). 
Because this evolution is sub-critical, the static configuration then completely disperses 
(frames 10-12) leaving flat spacetime behind. In a marginally super-critical case, a black 
hole containing essentially all of the mass of the static solution will form. 



non-linear dynamical system (with three continuous degrees of freedom) 
which exhibited many of the features seen in Type II critical collapse. 

6. Type I Behaviour: Einstein- Yang-Mills Collapse 

In contrast to the recently discovered Type II transitions, Type I black-hole 
transitions are, of course, very familiar to astrophysicists in the context of 
the instability which generically arises in sequences of static stellar models 
parametrized, for example, by central density. A somewhat different Type I 
transition has recently been observed in the spherically symmetric collapse 
of an SU(2) Yang-Mills field [35]. The equations of motion for this system 
are very similar to those for the EMKG model, with a single function, 
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W(r,t) (the "Yang-Mills potential") playing the role of the scalar field. 
However, largely due to the non-trivial vacuum structure of the theory, 
the phenomenology here is richer than in the scalar case. In particular, 
Bartnik and Mckinnon [36] showed that the Einstein- Yang-Mills (EYM) 
model admits a countable infinity of regular, static solutions, conveniently 
labeled by the number, n, of zeros of the static profile W n (r) (W must be 
±1 at the origin and at infinity). 

All of the Bartnik-Mckinnnon solutions, which may be viewed as result- 
ing from a balance between the attractive gravitational interaction and the 
repulsive YM self-interaction, were soon demonstrated to be unstable [37], 
and, in fact, within the specific ansatz in which they were originally con- 
structed, were shown to have precisely n unstable modes. This suggested to 
Bizoh that the n = 1 solution could be constructed as a Type I critical so- 
lution from the dynamical evolution of suitably constructed one-parameter 
families of initial data. This indeed turned out to be the case [35] , and Fig- 
ure 6 shows the evolution of a configuration near the Type I transition. As 
with the Type II solutions described above, the critical solution here has 
exactly one unstable mode in perturbation theory — tuning of initial data 
controls the effective level of that mode in the initial data, and as one ap- 
proaches the critical point, the unstable n = 1 Bartnik- Mckinnon solution 
persists for a proper time 

T <7ln|p-p*|, (27) 

where in complete analogy to the Type II analysis, the scaling exponent, 
a, is simply the reciprocal of the real part of the eigenvalue corresponding 
to the unstable mode. 

Interestingly, Type II behaviour is also seen in the EYM model [35] with 
A = 0.74 and 7 = 0.20, providing rather definitive evidence of the non- 
universality of mass-scaling exponents. Furthermore, considering suitable 
two-parameter families of solutions, S[p\,p2], one generically finds a co- 
existence point of the two types of threshold behaviour in the (pi,P2) plane. 
Phenomena analogous to all those described in this section have also been 
recently reported by Brady, Chambers and Goncalves [38] in the context of 
the collapse of a massive scalar field. 

7. An Irregular Critical Solution 

A final and novel example [39] of black-hole critical behaviour is provided by 
a family of solutions (family parameter e) of the EMKG system which has 
been independently discovered over the years by many authors [40, 41, 42], 
most recently by van Putten [43] . These solutions are spherically-symmetric 
and static, with an irregular origin, and have the feature that as the ad- 
justable parameter e is tuned to 0, the region of spacetime exterior to 
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Figure 7. Excitation of the growing mode of the static, irregular scalar field critical 
solution found by various authors [40,41,42]. and recently discussed by van Putten [43]. 
Here, the dynamical evolution of the quantity raa~ 1 d(j>/dt, which vanishes for the static 
solution, is plotted. The initial data contains a perturbation in the scalar field, which 
travels inward, bounces off the irregular origin between t = 18 and t = 28, then escapes 
to infinity. As it propagates through the region of the spacetime corresponding to the 
near-horizon region in a Schwarzschild spacetime, the perturbing pulse excites the sin- 
gle, growing mode of the static solution. Through judicious choice of perturbations, the 
growing mode can be excited with either sign — one sign leads to black hole formation, 
the other apparently leads to dispersal. Also note the blue-shifting and red-shifting of 
the perturbation (indicative of the strong-field nature of the static solution) as it travels 
inwards and outwards respectively. 



the effective support of the scalar field becomes a better and better ap- 
proximation of the Schwarzschild geometry. However, none of the solutions 
exhibit event horizons, and because of the difficulty in numerically treat- 
ing horizon-containing spacetimes, van Putten proposed that the solutions 
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might be useful for "mocking up" black holes. Unfortunately (from the 
point of view of that proposal), the solutions are unstable, and thus likely 
to be of little use as "approximate black holes" . However, as with the other 
threshold solutions described above, these solutions appear to be minimally 
unstable — that is, they seem to have a single unstable eigenmode in pertur- 
bation theory and, furthermore (see Figure 7), appear to sit at the threshold 
of black hole formation. 

8. Angular Momentum in Critical Collapse 

Very recent work by Gundlach [44, 45] addresses the important issue of 
the role of angular momentum at the threshold of black hole formation. 
In [44], non-spherical perturbations about the spherically-symmetric Evans- 
Coleman radiation- fluid (P = p/3) critical solution are considered, with the 
conclusion that there are no growing non-spherical modes. This strongly 
suggests that, at least for small departures from spherical symmetry, the 
Evans-Coleman solution will persist at the black hole threshold. In [45] 
Gundlach investigates the dominant non-spherical mode {I = 1), which 
contributes to the angular momentum L of the system and predicts that 
the overall amplitude of the spin should scale as 

Loclp-pT, (28) 

with [i ~ 0.80, but also that the spin vector will rotate in space as a 
function of p — p*, with a universal frequency lo ps 0.23. It will be extremely 
interesting to investigate these predictions, and related issues, using full 
simulations in axisymmetry. 

9. Quantum Mechanical Effects 

The nature of Type II critical solutions — particularly the fact that they de- 
scribe strong-field, gravitationally-mediated dynamics down to arbitrarily 
small scales — makes it natural to study the impact of matter-quantization 
on criticality. At the current time, the overall picture is murky, even given 
the restriction to semi-classical approaches. Here, there are at least two 
heuristic viewpoints. On the one hand, quantization will generally introduce 
a length scale to the otherwise scale-free dynamics at a Type II threshold. 
Thus, one might expect a mass gap to appear at threshold, and there is some 
evidence [46, 47, 48] for this. On the other hand, the predominant semi- 
classical effect which will be encountered is probably Hawking radiation, 
which, one could argue, should merely "renormalize" the classical competi- 
tion in the model, suggesting that the transition remains Type II [49, 50]. 
Despite the current confusion, one thing does seem clear — because Type 
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II solutions are naturally strong-field on arbitrarily small scales, they de- 
mand to be quantized in a way that the Schwarzschild solution, for example 
(where one must go in and, by hand, tune M to zero to achieve unbounded 
curvatures) does not. 



10. Conclusions 

In a few short years, the study of the threshold of black hole formation 
has developed into an active and (at least to some of us!) exciting sub- field 
of relativity. Although the initial results are now fairly well understood, 
it is clear that there remains an extremely rich phenomenology waiting 
to be explored, particularly in the context of non-spherically-symmetric 
collapse. The reader may rest assured that vigorous efforts to explore this 
new territory are currently underway. 
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